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Abstract 



An analytical method was developed to estimate errors in 
quantities depended on full one-port vector network analyzer 
(VNA) measurements using differentials and a complex differ- 
ential error region (DER) was defined. To evaluate the 
method, differences instead of differentials were placed over 
a DER which was then analyzed and compared with another com- 
monly used estimated error. Two real differential error in- 
tervals (DEIs) were defined by the greatest lower and least 
upper bounds of DER projections. To demonstrate the method, a 
typical device under test (DUT) was built and tested against 
frequency. Practically, a DER and its DEIs are solely based 
on manufacturer's data for standard loads and their uncer- 
tainties, measured values and their inaccuracies. 



Introduction 

In full one-port measure- 
ments with a VNA of real charac- 
teristic impedance Z G , a DUT 
with impedance Z has a reflec- 
tion coefficient p defined by 

p= (Z-Z 0 )/(Z + Z 0 ) 

and related to its measured 
value m by the bilinear trans- 
formation 

p = (m-D)/[M(m-D) + R] 

in terms of errors D, M and R 
[1] . This transformation can be 
uniquely determined from given 
distinct p n , n = 1, 2, 3 and re- 
spectively known m k , k = n [2] . 



Research 

We considered p n , m k as the 
elements of given ordered 
triples (A, B, C), (a, b, c), 

solved the resulting system 
and appropriately expressed 
its solution by 

F = X cC(B - A) 

D = X abC(A - B)/F 

M = X c(B - A)/F 

R= [l"l (A - B) (a - b) ]/F 2 

where X and |"l produce two 
more terms from the one 
shown, by rotation of the or- 
dered triple elements. These 
errors were then considered 
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as depended on the indepen- 
dent variables p n , m k . There- 
fore, their differentials 
were expressed in the same 
manner by 

dD = [n (a - b) X (B - C)BC dA 
+ X (b - c) 2 (B - A)(C - A) 

BC da] /F 2 

dM = [X (a - b)(c - a)(B - C) 2 dA 

- n (A - B) X (b - c)da]/F 2 
dR = {X [F + 2(a - b)B(A - C)] 

[(B - C) 2 dA n (a - b) 

- (b - c) 2 da n (A - B) ] }/F 3 

After that, the differential 
of p was expressed by 

dp = [-RdD - (m - D) 2 dM - (m - D)dR 
+ Rdm]/[M(m - D) + R] 2 

and was considered depended, 
through dD, dM and dR, on L = 
7 independent variables and 
their independent differen- 
tials: p n , n = 1, 2, 3 and m k , 
k = n or k = 0 with m 0 = m. 

The developed expressions 
were mechanically verified 
using a developed software 
program for symbolic computa- 
tions . 

Manufacturer's data for 
standard loads used in full- 
one port VNA measurements are 
substituted in p n , and for 
their uncertainties in dp n . 
Since Z 0 is real, the domain 
of each p n is the closed unit 
circle [3]. For |p n | = 0 or 1, 
care must be exercised to re- 
strict its differential value 
onto its domain. The VNA mea- 



surements have specified 
bounded ranges for their mo- 
dulus and argument, so that 
the domain of each m k is a 
bounded circular annular with 
its center at the origin O of 
the complex plane. Measure- 
ment data are substituted in 
m k and manufacturer's data 
for measurement inaccuracy in 
dm k . Uncertainty and inaccu- 
racy data outline domains for 

dp n and dm k . If z = | r | e- 1 ^, stands 
for any of the independent 
variables and dz for its dif- 
ferential then the contribu- 
tion of dz to dp is a summa- 
tion term of the form Wdz, 

with W = | U | e^ v , so that 
Wdz = |U|e j(v + <p) d|r| 

+ |u|e j(v + (p+ TT/2) | r | dtp 

where W is in fact a known 
value of the respective par- 
tial derivative and d|r|, dtp 
are the independent real dif- 
ferentials of the complex dz 
in polar form. Each expres- 
sion Wdz outlines a contour 
for a partial DER around O. 
If z ^ 0, the partial DER is a 
parallelogram with perpendic- 
ular sides d | r | and | r | dcp, 
stretched or contracted by |U| 
and rotated by (V + tp) around 
O. If z = p n = 0, the partial 
DER is a circle with radius 
| U | d | r | . Accordingly, a DER 
is the sum of either L paral- 
lelograms or ( L - 1 ) parallel- 
ograms and 1 circle. DER is 
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then a convex set with con- 
tour either a polygonal line 
with 4L vertices at most, or 
a piecewise curve composed of 
4( L - 1) line segments and 4(L 
- 1) circular arcs at most. 
The greatest lower and least 
upper differential error 
bounds are the end-points of 
DEIs for the real and imagi- 
nary parts of dp and result 
from the projections of DER 
for p on the coordinate axes. 
These conclusions can be ge- 
neralized for any other quan- 
tity directly or indirectly 
depended on all, some or just 
one of the above independent 
variables and their differen- 
tials. Thus, the quantity has 
an L-term DER, where 7>L>1. 
For example, the impedance Z 
of a DUT has the 7-term DER: 

dZ = 2Z 0 dp/(l - p) 2 

Results 

All of the following data 
are specified by manufactur- 
ers of the parts for our mea- 
surement system. This system 
operates from 1 to 1300 MHz 
with 100 Hz PLL stability and 
consists of a type-N Z 0 = 50 Q 
network analyzer, a number of 
support instruments and a set 
of standard loads. The stan- 
dards are: a short circuit A, 
a matching load B with re- 
flection coefficient 0.029 
and an open circuit C with 
reflection coefficient 0.99 
and phase accuracy ±2°. In 
the absence of manufacturer's 



data for A we considered its 
uncertainty equal to that of 
C. So, the following values 
were substituted in the de- 
veloped expressions: 

A = -1, 0 < d | A | <0.01, -180° < 
dtp A < -178° or 178° < dtp A < 180°, 

B = 0, | dB | = 0.029, 

C = 1, -0.01 < d | C | <0, -2° < dtpc 
< + 2 ° 

The annular domain for m k 
of VNA is specified from 0 to 
-70 db in modulus and ±180 
degrees in argument. Measure- 
ments m k result with a deci- 
mal floating point precision 
of 4 digits, for both modulus 
and argument. We consider the 
modulus and argument of dm k 
equal to ±1/2 of the unit in 
the last place of the corre- 
sponding mantissa in modulus 
and argument of m k Conse- 
quently, our system produces 
a DER, either for p or Z, as 
a sum of ( L - 1 ) = 6 parallelo- 
grams and 1 circle, with a 
contour of (4L + 4L) = 48 ver- 
tices at most. 

A suite of developed soft- 
ware applications: (i) con- 

trols the system and collects 
the data in terms of fre- 
quency using the IEEE-488 
protocol, (ii) processes the 
collected data and computes 
the vertices of DER and the 
end-points of its DEIs (iii) 
sketches pictures of DER for 
p and Z in terms of the fre- 
quency steps and make a film 
using them as frames. 



THURSDAY 31 OCTOBER 2013 



V4-19 



FUNKTECHNIKPLUS # JOURNAL 




A/. I. YANNOPOULOU, P.E. ZIMOURTOPOULOS 
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Input Resistance R, Ohm 



Fig.l: A typical DER for the impedance Z 




0.2 

Fig. 2: DER for the reflection coefficient p and for its 
associated impedance Z against frequency 
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Fig. 3: Greatest lower and least upper differential error 
bounds for resistance R and reactance X against frequency 



A typical resistor with a 
nominal DC impedance of 50 Q 
±20% tolerance was soldered 
on a type-N base connector 
and enclosed in an aluminium 
box to serve as a simple DUT 
for testing its Z from 2 to 
1289 MHz in 13 MHz steps. The 
center frequency f c = 639 MHz 
was chosen to reveal the de- 
tails of the proposed method 
in Fig.l, where the contour 
of a typical DER for Z is 
outlined with small circles 
as its vertices. This contour 
surrounds that of the 4-terms 
DER due to inaccuracy of mea- 
surements (1) and that of 3- 



terms DER for the uncertainty 
of loads (2). A properly cir- 
cumscribed rectangle of DER 
shows graphically how the 
DEIs for R and X result. The 
commonly used error from the 
matching load only is shown 
as a dotted circle. This is 
in fact a 1-term DER which is 
surrounded from the contour 
of the DER by a factor of 
about 125% to 185% in all di- 
rections. Finally, in the 
same figure, 2 7x2 differences 
AZ resulting from the same 
dp n and dm k , dense enough to 
appear as stripes, were 
placed over DER to compare 
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them with differential dZ 
values. Notably, almost all 
of AZ values are belong to 
DER while the computation 
time for these AZ exceeds 
that for DER by more than one 
order of magnitude. To demon- 
strate the method, a set of 
selected DER frames for p and 
Z are shown in Fig. 2, as 
beads on space curved fila- 
ments against frequency. 



Finally, the computed DEIs 
for R and X are shown in 
Fig. 3 against frequency. 

Conclusion 

The proposed method may be 
efficiently used in the same 
way, to successfully estimate 
errors in any quantity de- 
pended on full one-port vec- 
tor network analyzer measure- 
ments . 
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